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and most were confirmed as probable high-redshift quasars by supplementing the SDSS 
data with J and K measurements. The quasars possess strong, broad Lya emission 
lines, with the characteristic sharp cutoff on the blue side produced by Lya forest ab- 
sorption. Three quasars contain strong, broad absorption features, and one of them 
exhibits very strong Nv emission. The amount of absorption produced by the Lya 
forest increases toward higher redshift, and that in the z = 5.27 object (Da ~ 0.7) is 
consistent with a smooth extrapolation of the absorption seen in lower redshift quasars. 
The high luminosity of these objects relative to most other known objects at z }t 5 makes 
them potentially valuable as probes of early quasar properties and of the intervening 
intergalactic medium. 



Subject headings: quasars: general - surveys 



1. Introduction 



Since the identification of the first quasar redshift (3C 273, Schmidt 1963), quasars have been 
at the forefront of modern cosmology. With luminosities tens or hundreds of times higher than 
those of galaxies, quasars are a powerful probe of the distant primordial universe. Over the past 
decade approximately two hundred objects above redshift four have been discovered, and there is a 
growing consensus that the number density of luminous quasars peaks between redshifts of two and 
three and steeply declines out to the limits of current measurements (z « 4.5; see Warren, Hewett, 



& Osmer 1994; Schmidt, Schneider, fc Gunn 1995 ; Kennefick, Djorgovski &; de Carvalho 1995 ). 



Recent studies ( Fan et al. 1999| ,2000) have dramatically increased the number of known quasars 
with redshifts larger than 4.5, opening the possibility of investigating the quasar luminosity function 
at redshift five and beyond; this information will determine whether the number density of quasars 
continues to decline with increasing redshift, or whether theoretical models (e.g., Haiman & Loeb 



1998) that predict a significant number of quasars at z > 5 are correct. 



The majority of high-redshift quasars have been identified by optical color selection. As the 
result of intergalactic absorption, the flux in the spectral region shortward of Lya in distant objects 
is significantly attenuated. The onset of the Lya forest and the Lyman break can be detected with 
broad-band imaging. Multicolor optical surveys ( Warren, Hewett, &: Osmer 1991 , Irwin. McMahon 



and Hazard 1991, Djorgovski et al. 1999) have proven effective in identifying z > 4 quasars, but 
as most such surveys lack information redward of the I band, they have difficulty distinguishing 
between quasars at redshifts larger than ~ 4.8 and cool stars. Recently, the commissioning data of 
the Sloan Digital Sky Survey (SDSS, York et al. 2000|) have led to the identification of new high- 
redshift quasars at an unprecedented rate: approximately 45 quasars at z > 3.6 (two at z ~ 5, Fan 



et al. 199S ,2000; [Schneider et al. 2000 ), have been published in the past two years. These quasars, 
at Mb ~ —27, are at the luminous end of the quasar luminosity distribution. (Throughout this 
paper cosmological properties are calculated assuming H Q = 50 km s" 1 Mpc -1 and q = 0.5.) 
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In addition to wide area surveys, color selection for high-redshift objects has been applied to 



small, deep fields. Three galaxies at z > 5.6 have been spectroscopically confirmed (Hu, Cowie, fc 



McMahon 1998| ; |Weymann et al. 1998| ; |Hu, McMahon, fc Cowie 1999| ), and recently [Stern et al. 



(2000)| found an AGN (M B ~ -22.7) at z = 5.5. The Hubble Deep Field has led to the discovery of 



high-redshift candidates (z ~ 6-10, [Lanzetta et al. 1999 ; Chen et al. 1999| ); these objects, however, 



are orders of magnitude less luminous than quasars, and their spectral properties are difficult to 
study even with the world's largest telescopes. 

Although the z = 5 barrier was broken more than two years ago, objects at these extreme 
redshifts are still rare enough that each new example is a potentially valuable probe of the very 
early universe. This is especially true in the case of the most luminous quasars, since follow-up 
spectroscopy of these relatively bright objects can reveal information about the physical conditions 
in early quasars and about the state of the intergalactic medium at very high redshifts. The SDSS 
collaboration has undertaken extensive efforts to search for high-redshift quasars. In this paper, 
we report the discovery of five quasars at z > 4.6, of which the most distant is at z = 5.27. 



2. Selection of Quasar Candidates 



The SDSS (York et al. 2000) utilizes a wide-field camera with 54 CCDs QGunn et al.T998p 



mounted on a dedicated 2.5m telescope at the Apache Point Observatory (APO), New Mexico, to 
survey ~ tt steradians of the sky around the Northern Galactic Cap. CCD images in five broad 



optical bands (vf, g' , r', i' , z', centered at 3540 A, 4770 A, 6230 A, 7630 A and 9130 A ; [Fukugita et 



al. 1996 ) yield a nominal 5a detection of point sources in AB magnitudes of 22.3, 23.3, 23.1, 22.3, 
and 20.8, respectively. The commissioning data have so far covered ~ 600 square degrees, mostly 
near the equatorial region, \S\ < 1.5°. 

With five bands and a spectral resolution of R ~ 3, the SDSS imaging data can distinguish 
quasars from stars in a broad redshift range ( [Fan 1999|) . At 3.5 < z < 5, quasars lie well away from 



the stellar locus in the g'r'i'z' color space, due to the large equivalent width of the Lya emission line 
and the significant absorption produced by the Lya forest and Lyman limit systems. At redshifts 
between 4.4 and 5, the r' — i' color becomes large, and i 1 — z' is near zero. At z ^ 5, the flux 
drop shortward of the redshifted Lya emission line affects the i'-band magnitude, and the i! — z' 
color increases with redshift. At this point the quasar track quickly approaches the red end of 
the stellar locus in the r'i'z' color diagram; the rise in contamination by very cool stars requires 
that additional discriminators be added to aid the quasar selection. For quasars, the underlying 



continuum longward of Lya can be approximated with a power law ([Francis et al. 1991| ; [Schneider, 



Schmidt fc Gunn 1991| ; |Zheng et al. 19971 ) of f v ecu , leading to a small color differences in 



bands located redward of the Lya emission line (i.e., blue band — red band ~ 0). This small color 
difference contrasts sharply with the colors of most cool stars ( Leggett et al. 1996j , 2000), whose 
flux rises rapidly towards longer wavelengths. 
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Our target selection is based on three regions in SDSS color space: (1) r* — i* > 1.35 and 
i* - z* < 0.3; (2) r* - i* > 2 and i* - z* < 0.7; or (3) z'-band detection only, i.e. z* < 20.8 and 
the detection in the other four bands is below 5a. If an object's r' or i' magnitude is below the 
respective 5a detection level, we use the latter in calculating the color in case (1) and (2). (Note 
the * superscript for the magnitudes; the measurements reported here are based on a preliminary 
calibrations of the SDSS photometric system.) In addition, an object must be classified as a point 
source by the SDSS processing software and have z* < 20.8 to be included as a quasar candidate. 
We applied these selection criteria to ~ 200 square degrees of SDSS imaging data acquired in 1999 
March and 2000 February. 

The selection criteria used here differ from those employed by Fan et al. (1999, 2000) to identify 
z > 4 quasars. Whereas Fan et al. required i* < 20.2, candidates in this paper can be undetected 
in the i' band. The flux measurements for this paper's candidates often have lower signal-to- 
noise ratio, and a sample drawn from these selection criteria will naturally be more susceptible to 
contamination from non-quasars than those in Fan et al. The color cuts in this paper are closer 
to the stellar locus than the Fan et al. criteria; this will also resultin a reduction in the selection 
efficiency. 

We obtained IR photometry of quasar candidates in the J and K bands at the NASA Infrared 
Telescope Facility (IRTF) at Mauna Kea, Hawaii. The observations were taken on 2000 March 9-12 
using the NSFCam equipped with a 256x256 InSb array. The plate scale was / .'30 pixel -1 , and the 
seeing was ~1".2 in the if -band. Each selected target was imaged with standard dithering technique 
with total exposure times of 7 min and 5 min in J and K, respectively. Images of IR standard stars 
were taken throughout the observations to monitor the magnitude zero points. Typical photometric 



errors are in the 0.05-0.1 magnitude range, depending on the brightness of objects; see Tsvetanov 



et al. (2000) for details of the observations and data calibration. We observed ~ 20 known SDSS 



quasars with redshifts greater than 4 to empirically calibrate our selection technique. Our additional 
color constraints are z* — J < 1.5 and J — K < 1.8. Approximately ten such objects whose colors 
closely resembled those of known quasars were selected. Note that four of the five quasars are at 
redshift smaller than 5, and they can be selected by the SDSS data alone ([Fan et al. 1999| , 2000). 



3. Spectroscopic Observations 

Spectroscopic follow-up observations of SDSS high-redshift quasar candidates were carried 
out in 2000 February and March, with the Digital Imaging Spectrograph (DIS) of the APO 3.5m 
telescope. The DIS is a double spectrograph; for high-redshift quasars, only the red part of the low- 
resolution spectrum, covering the wavelength range 5400 A to 10000 A at 13 A resolution, contained 
any useful signal. The exposure time for each object was 30 minutes; even with a limited signal-to- 
noise ratio in the spectra, the redshift identifications, based on the strong, asymmetric Lya emission 
line and absorption produced by the Lya forest, are unambiguous. The observations in February 
were made before the IRTF run. Of the six candidates, only one, SDSS 1129—0142, turned out to 
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be a quasar. The J2000 coordinates are given in the object name (format hhmmss.ss+ddmmss.s; 
see Fan et al. 1999| .) For brevity, we have shortened the names to SDSS hhmm+ddmm throughout 



the text.) In late March we observed five SDSS/IR candidates, and two, SDSS 1021—0309 and 
SDSS 1208+0010, are quasars. 

After the initial identification with the APO data, additional spectroscopic observations of the 
three quasars were obtained in 2000 April with the Low Resolution Spectrograph (LRS; pill et aT 



2000[ [Schneider et al. 2000| ) at the prime focus of the Hobby-Eberly Telescope ( |Ramsey et al. 1998[) 
at McDonald Observatory. The LRS configuration of a 2" slit, 300 line mm -1 grating, and OG515 
blocking filter produced spectra with a wavelength coverage of 5150 A to 10,150 A at a resolution 
of 20 A. The exposure times were typically 30 minutes. 

Two additional SDSS/IR quasar candidates were observed with the Keck 10m telescope on 
April 5-6, and both of them, SDSS 1451—0104 and SDSS 1122—0229, are quasars. The spectra 
were taken with the Echellette spectrograph and imager (ESI, Epps Sz Miller 1998j ) on the Keck 



Observatory 10-m telescope. The ESI was used in high dispersion mode which covers the wavelength 
range of 4000 to 11000 A with 11 kms -1 resolution. The quasars were viewed through a 1" wide 
slit oriented at the parallactic angle. The exposure time is 20 minutes each. The z = 5.27 quasar, 
SDSS 1208+0010, was also observed for 30 minutes. The spectra were flux calibrated relative to the 
standard star G191B2B. The spectra were extracted and reduced using standard IRAF programs, 
and binned to 3. 85 A. 

Only one of the five "non-infrared" candidates, observed in February, is a quasar. Of the seven 
SDSS/IR candidates observed, four are quasars, while the other are late-type stars. This result 
tentatively suggests that IR selection may significantly improve the selection efficiency, but clearly 
a larger sample is needed to confirm this conclusion. Table 1 lists the optical/IR photometric 
measurements. 

The quasar spectra are displayed in Fig. 2, with prominent spectral features marked. All the 
spectra have been placed on an absolute flux scale by matching the i* magnitudes in Table 1 with 
the i* magnitudes synthesized from the spectrum. The Keck and HET spectra of SDSS 1208+0010 
reveal the Civ emission at ~ 9500 A, which is not clear in the data taken at APO. 



Discussion 



Table 2 contains the redshift (see notes below for the specifics of the measurements for each 
object), AB1450 (the AB magnitude of the quasar at 1450 A in the rest frame, corrected for Galactic 
absorption), the power law index of the continuum, the continuum depression due to absorption 



in the Lya and Ly/3 regions (Da and Db; see pke &; Korycansky 1982j ), and the absolute B 



magnitude of the quasars. The Galactic extinction is calculated using the reddening map of Schlegel 



Finkbeiner Davis (1998) . The luminosities for the quasars were calculated assuming that the 
continuum power law slope from the far ultraviolet to the optical was —0.5 (Schneider, Schmidt 
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& Gunn 1991). The individual power law slopes, and hence the depression estimates, are quite 
uncertain given the limited baseline available in the spectra and the presence of BAL features. All 
are moderately luminous quasars, and the spectra of three contain BAL features. Each spectrum 
contains at least one excellent candidate for a damped Lya system, and there is no detected flux 
below the rest frame Lyman limit in any of the quasars. 

As shown in Figure 2, the extent of intergalactic absorption shortward of the redshifted Lya 
emission increases from z = 4.67 to z = 5.27. This is reflected in the continuum depression values in 
Table 2. However, the residual flux can be seen in all the spectra, particularly around the Ovi+Ly/3 
feature, suggesting that the intergalactic medium is not completely opaque. Our measurements are 
consistent with the known distributions of the Lya forest absorption as derived at lower redshifts 
(press fc Rybicki 1993|) . 

The emission features in these quasars are common among AGN at lower redshift. Three of 
the the quasars exhibit significant BAL features. Finding charts for the five quasars are given in 
Fig. 1. 

Notes on the quasars: 

SDSS 1021—0309 {z = 4.70): The sharp split of the Lya and Nv lines may suggest unusually 
strong Nv emission. A simple fit to a 300A- region longward of the Lya cutoff yields a Nv/Lya 
ratio of 0.64 for the narrow features. If the split is a result of Nv absorption, the corresponding 
wavelength for the C iv absorption counterpart should be centered at ~ 8710A. While a significant 
C iv BAL trough is present in the spectrum, it is centered at ~ 8560A. The redshift determined by 
the Ol, SilV+Oiv], and Civ emission lines are consistent with each other at the 0.005 level. The 
Lya and Nv emission shows redshifts that are consistent (with slightly larger errors) with the other 
lines. 

SDSS 1122—0229 (z = 4.80): This object has the strongest emission lines among the five. The 
measurements of C iv, Ol, Nv and Lya emission yields a redshift of z = 4.80 ± 0.03. The Lya 
absorption cutoff is not as sharp as the others, stretching ~ 60A with several narrow absorption 
troughs. 

SDSS 1129—0142 (z = 4.85): This quasar is slightly more luminous than 3C 273, which in the 
adopted cosmology is Mb = —27.0. It displays a number of spectacular BAL features; because of 
this, many of the properties given in Table 2 contain large uncertainties. Strong, broad troughs of 
C iv and Siiv+Oiv] dominate the spectrum, and there is a suggestion of the presence of Ol and 
Nv absorption features. The redshift is based on assigning a rest wavelength of 1219 A to the 
peak of the Lya emission line (see [Schneider, Schmidt Gunn 19911 ). The observed values of the 
continuum depression are exceptionally large; this almost certainly arises from significant intrinsic 
absorption. 

SDSS 1208+0010 (z = 5.27): Very strong, relatively narrow Lya and Nv emission dominate 
the spectrum; this spectrum bears an uncanny resemblance to that of the z = 4.04 quasar 
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PC 0910+5625 (ISchmidt, Schneider, fe Gunn 1987] ). The Nv,Oi, and Civ lines yield a consis- 



tent redshift; the peak of the Lya feature occurs at 1218 A, typical for quasars at redshifts above 
four ( [Schneider, Schmidt fc Gunn 1991 ). The depression due to the Lya forest (see Table 2) is 



quite large, but not unusual for this redshift, suggesting that there is not a dramatic change in the 
character in the Lya forest from lower redshifts to at least z ~ 5.2. 

SDSS 1451-0104 (z = 4.67): The redshifts determined from the peaks of Civ and Siiv+Oiv] 
match that of the Lya edge within 0.005. A significant, broad absorption trough is present between 
~ 8250 - 8550A. 

The discovery of these quasars once again demonstrates the ability of the SDSS to effectively 
identify z > 4.6 quasars, and extends the SDSS redshift range to well beyond five. By supplementing 
the SDSS measurements with J and K photometry, we have been able to efficiently identify (success 
rate of 50%) high-redshift quasars in magnitude/color space regions that are fainter and closer 
to the stellar locus than were presented in Fan et al. (1999,2000); it is likely that IR photometry 
will be a valuable tool in the search for faint, z > 5 quasars. Our IR photometry is only a test, and 
the results do not constitute a complete sample. To date the SDSS has imaged but a few percent 
of the planned survey area; based on the results to date, the complete survey should contain well 
over a hundred z > 4.7 quasars found with well-defined selection criteria. 

The Sloan Digital Sky Survey (SDSS) is a joint project of The University of Chicago, Fermilab, 
the Institute for Advanced Study, the Japan Participation Group, The Johns Hopkins University, 
the Max-Plan ck-Institut fur Astronomie, Princeton University, the United States Naval Observa- 
tory, and the University of Washington. Apache Point Observatory, site of the SDSS, is operated 
by the Astrophysical Research Consortium. Funding for the project has been provided by the 
Alfred P. Sloan Foundation, the SDSS member institutions, the National Aeronautics and Space 
Administration, the National Science Foundation, the U.S. Department of Energy, and Monbusho, 
Japan. The SDSS Web site is http:/ / www.sdss.org/] . 
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Fig. 1. — Finding charts of the five quasars. All the charts are i'-band images from the SDSS. 
North is up, East to the left, and the individual charts are 240" on a side. 

Fig. 2. — Spectra of the five quasars. The wavelengths of identified and expected major emission 
features are marked. The telescope names for individual spectra are marked at the upper-left corner 
of each panel. The fluxes have been normalized to the respective SDSS i*-band magnitudes. The 
spectral resolution is 13A for the APO data, 20 A for the HET data, and Keck spectra have been 
rebinned to 3.85 A. 
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Table 1. Optical Positions, SDSS a and Infrared b Photometry 



Object c 


u* 


9* 


r* 


i* 


z* 


J 


K 


E(B-V) 


SDSSp J102119.16-030937.2 


23.72 


25.82 


21.76 


20.09 


20.02 


18.77 


17.08 


0.042 




±0.55 


±0.47 


±0.09 


±0.03 


±0.10 


±0.07 


±0.05 




SDSSp J112242. 98-022905.1 


23.73 


24.77 


22.22 


20.38 


20.47 


> 19.5 




0.055 




±0.57 


±0.50 


±0.12 


±0.04 


±0.15 


±0.10 






SDSSp J112956.10-014212.4 


24.06 


25.26 


22.02 


19.64 


19.51 


17.51 


16.09 


0.072 




±0.55 


±0.37 


±0.10 


±0.03 


±0.07 


±0.05 


±0.05 




SDSSp J120823.82+001027.7 


24.39 


24.77 


22.75 


20.79 


20.72 


19.43 


18.10 


0.024 




±0.65 


±0.31 


±0.08 


±0.27 


±0.35 


±0.10 


±0.10 




SDSSp J145118.77-010446.2 


24.40 


24.45 


22.65 


20.70 


20.53 


19.43 


18.17 


0.044 




±0.45 


±0.47 


±0.16 


±0.04 


±0.13 


±0.10 


±0.10 





and 2* bands are approximately 24.00, 24.85, 24.55, 24.20, and 



a Preliminary SDSS calibration; values are expressed as asinh magnitudes (Lupton, Gunn, 
asinh magnitudes for zero flux in the u* , g* , r* 
22.40, respectively. 

b In UKIRT (Vega-based) magnitude system. 



Szalay 199S). The 



c The J2000 coordinates are given in the object name (hhmmss.ss±ddmmss.s). 



Table 2. Quasar Properties 



Object 


z 


AB1450 


a a 


D A 


D B 


M B h 


SDSSp J102119.16-030937.2 


4.696 ± 0.004 


20.25 


-1.6 


0.58 


0.76 


-26.3 


SDSSp J112242. 98-022905.1 


4.795 ± 0.004 


20.58 


-1.2 


0.64 


0.79 


-26.0 


SDSSp J112956.10-014212.4 


4.85 ± 0.03 


19.22 


-1.3 


0.84 


0.97 


-27.4 


SDSSp J120823.82+001027.7 


5.273 ± 0.004 


20.47 


-0.7 


0.71 


0.81 


-26.3 


SDSSp J145118.77-010446.2 


4.672 ± 0.004 


20.42 


-0.7 


0.71 


0.86 


-26.2 



a Spectral energy index (/„ oc v a ); typical uncertainties are several tenths. 

b Calculated assuming Ho = 50, qo = 0.5, and a spectral energy index between the ultraviolet and blue band of 
-0.5. 
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Fig. 2.- 



